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O ■ ABSTRACT 

CD : 

Q . The dependence of the number density of the bright QSOs at different redshifts 

{nQso{.z)) on initial power spectrum is studied. It is assumed that QSO phenomenon 
is an early short term stage of evolution of massive galaxies with M > 2 x lO^^/i^^M©. 
The duration of such QSO stage which is passed through by fraction a of galaxies is 
^ . determined by means of minimization of the divergence of the theoretical number density 

lO \ of QSOs at different redshifts for specified initial spectrum from observable one [0. It 

is shown that the nearest number densities of QSOs at 0.7 < z < 3.5 to observable 
C^ ■ ones are obtained for the tilted CDM model (fib = 0.1, ra = 0.7). The QSO stage lasts 

O , ~ 7 X 10^/a years and begins soon after the moment of rise of the first counterflow in 

Q^ \ collisionless component and shock wave in gas. 

The possibility of the reconstruction of initial power spectrum on small scale on the 
Q^' base of the observable data on number density of QSOs at different z is considered too. 

Q . Such reconstructed spectrum in comparison with standard CDM has steep reducing of 

^ \ power at fc > 0.5hMpc~^. 



in 



^ 



1 



Subject headings: cosmology: initial power spectrum - dark matter - quasars 



1. Introduction 

The problem of formation of galaxies and large 
scale structure of the Universe has an astrophysi- 
cal aspect consisting in search of connection between 
properties of luminous objects and dark matter halo, 
as well as cosmological aspect consisting in determi- 
nation of the initial power spectrum of density fluc- 
tuations. Obtaining such spectrum from first princi- 
pals would essentially simplify the solving of prob- 
lem on the whole. Unfortunately, today we can 
not surely say which the postinflation spectrum is: 
scale invariant Ppi = Ak" with n = 1, tilted with 
n ^ 1, or more complicated. The nature of dark mat- 
ter (baryonic, coUisionless cold dark matter, futher 
CDM, hot dark matter, futher HDM, or their mix- 
ture H-I-CDM, or other), and values of cosmological 
parameters h = -, „„,. J^/° .,, „^ {Hq- Hubble constant), 
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- mean density of baryons) and cosmological constant 
A remain uncertain. More preferable from theoretical 
point of view spectra (standard CDM model, tilted 
CDM, CDM-f A or hybrid H-I-CDM) give predictions 
consistent with the observable data marginally only 
||^|[. Futhermore, the different spectra normalised 
to the COBE r.m.s. cosmic microwave background 
(CMB) anisotropy predict close characteristics of the 
large scale structure of the Universe, e.g. AT/T at de- 
gree angular scale, large scale peculiar velocity field 
etc. [|[ |[ ^ so that it is impossible to distinguish 
them. Therefore, it is very important to build up 
the new cosmological tests based on extension of a 
class of objects with predicted observational manifes- 
tations sensitive to initial power spectrum on small 
scale. Here we shall investigate the dependence of 
the number density of quasars (QSOs) at different 
redshifts against initial power spectra. To make it 
possible we suppose that a QSO is the early short 
term stage of evolution of a massive galaxy and ap- 
pears in the corresponding peaks of initial random 
Gaussian field of density fluctuations. The number 
density of already collapsed peaks of certain scale at 
different redshifts z strongly depends on an amplitude 
and slope of the power spectrum of such field at small 
scales (fc - 0.1 - lOhMpc-^). 

Such approach was applied independently by some 
authors. For example, Cen et al. |y have calcu- 
lated the comoving number density of nonlinear ob- 
jects as a function of redshift z using Press-Schechter 
approximation [lO[. Their results show that observ- 



able number density of quasar-like objects between 
2 < z < 4 can be explained by tilted CDM with 
n > 0.7 and mass of host galaxies < 5 10^^ Mq. 
Blanchard et al. ]ll| tested HDM models and have 
shown the disability of these models to explain the 
observable concentration of QSOs at z ~ 3 — 4 when 
primordial scale-invariant spectrum (n = 1) is nor- 
malised to the COBE r.m.s. cosmic microwave back- 
ground anisotropy |12|. Kashlisky [Q tested the hy- 



brid H-I-CDM and modified CDM models by num- 
ber density of distant QSOs at z = 4.5. The num- 
ber density of the peaks-precursors of massive galax- 
ies {M > 10^^h~^MQ), which passed a quasars stage 
soon after their central regions had collapsed, was cal- 
culated for the different spectra. The mean number 
density of QSOs predicted by modified CDM model 
with nh = 0.1, 0.2, 0.3 and H-hCDM was found to 
be less than observable one Uq^q « iOOh^Gpc~^ at 
z = 4.5. Nusser & Silk ||lj] for finding nqsoi^) in 
redshift range 2-5 for given spectrum used the Gaus- 
sian peaks formalism by Bardeen et al. |Iq| and 
shown that this formalism is more preferable than 
Press-Schechter one. They assumed the number den- 
sity of very bright QSOs {L > lO'^'^erg/s) to be 
equal to the number density of peaks with masses 
M > 2 X 10^^h~^MQ within certain height interval 
corresponding to the physically substantiated dura- 
tion of quasar stage (« 4 x 10^ years). They have 
concluded that tilted CDM model with n — 0.8, 
H+CDM {ncDM ~ 0.7, nHDM ~ 0.3) and scaled 
CDM with 0.2 < O/i < 0.3 are compatible with the 
mean abundance of quasars at 2 < z < 5. Ma & 
Bertchinger pq | have carried out the numerical sim- 
ulation of formation of elements of large scale struc- 
ture of the Universe with masses 10^^ — IO^^Mq in 
hybrid H-I-CDM models and analysed the possibil- 
ity of explanation of the observable number density 
of QSOs at 2.2 < z < 4.5 in the framework of the 
Press-Schechter formalism. They have shown that 
H-I-CDM model with ^hdm = 0.2 predicts nqsoiz) 
which matches the observable one if quasars are cre- 
ated soon after collapse of nearly spherical peaks. 

The goal of the paper is an improvement of this 
approach to the explanation of the observable data 
on abundance of the bright QSOs in wider range 
z — 0. 7-^4. 7 ||l|, |l7|, testing of the preferable ini- 
tial spectra, which marginally match other observ- 
able data on large scale structure of the Universe, 
search of spectrum momenta to reproduce the observ- 
able data on nqsoiz) and construction the suitable 



phcnomcnological spectrum on small scale. 

2. The main assumptions and method 

We explore the distribution of quasar abundance 
over redshift nqsoiz) in the framework of the the- 
ory of formation of the large scale structure of the 
Universe in the peaks of random Gaussian field of 
scalar density fluctuations as consequence of gravita- 
tional instability on the Einstein-de Sitter cosmolog- 
ical background. In this approach galaxies and ob- 
jects of other scales are formed in the peaks of den- 
sity fluctuations which have an amplitude S = — > 
I'thO'oiRf), where (Tq - r.m.s. amplitude of the fluctua- 
tions of certain scale Rf, vth - some threshold height 
which corresponds to minimal amplitude when ob- 
jects of this scale forms still. 

In the theory of random Gaussian fields ||l^, Q 
a comoving number density of the peaks with height 
v > vth is 



n{vth) 



p{v; Vth, q) Npk{v) dv, (1) 



where p{v]Vth-,q) = TMvJvyi' threshold function, 
^pk(y) is differential number density of peaks in the 
amplitude range (y, v + dv), which depends on v and 
spectrum momenta ctoi fi: f2 (formulae 4.3-4.5 and 
A15 in the paper by Bardeen et al. |15[ ). 

We suppose that certain fraction of massive galax- 
ies a\ passes through quasar stage of evolution. Maybe, 
only some fraction of them ai is visible by Earth 
observer. Therefore, the fraction of galaxies passing 
through quasars stage and visible by Earth observer 
as QSOs is a = aia^ < 1. If we adopt the black hole 
model of quasars, then the minimum mass associated 

10" 



IQI^Mq 1 14 



H 



with a host galaxy is 
For selection of galactic scale peaks with mass M > 
2x 10^^/i~^Mo we smoothed the initial random Gaus- 
sian density field (5(lzf ) by the Gaussian filter function 

- ^{lt;Rf) = ^,j2y. lS{P) e ''"/ d^7 with 
Rf — 0.35ft.^"'^Mpc. It means that the power spec- 
trum of the smoothed field is P(fc; Rf) = P(fc)e"'''^/ , 
where P{k) =< SkS'^ > is initial power spectrum, 
dk is Fourier transformation of the 6(lif): 6(jf) = 



■fS.e 



^ d"^ k . Those among them which 
> 



(27i-)i-5 

have an amplitude 5{lif,Rf) 



"th'^o, 



where CTq 



(2^ 



l^Jdkk^P{k) 



-k^Ri 



is r.m.s. amplitude of den- 



sity fluctuations, evolve to massive galaxies. Their 
concentration is known from CfA cataloque |2^ 



(4.6/i ^Mpc) ^. Using eq.(l) for massive galaxies 



'^i^th^'f) =n 



obs 



(2) 



one parameter of threshold function (e.g. i^f^) can be 
found. 

As follows from Tolmen model, a homogeneous 
spherical-symmetrical adiabatic fluctuation on the 
Einstein-de Sitter cosmological background collapses 
at 

'to 

tc 



1 = Q.59aoVc - 1, 



(3) 



where to — 1-3 x 10^" years is present cosmological 
time (here and further we assume Hq = 50 ^^^ ), tc 
- the moment of collapse, Vc - corresponding height of 
peak. At this moment tc appearence of the first coun- 
terflows in dark matter and generation of shock wave 
in baryon component take place. Massive black hole 
in central region of such peak may be formed during 
^ 10^ years as it follows from the physical models 
of quasar mechanism (see, for example, [pO|). There- 
fore, we suppose that QSO stage of galaxy evolution 
begins soon after collapse of central region of peak i^ 
or later by certain time interval At called a time de- 
lay. Duration of the QSO stage or quasar lifetime is 
tqso- 

In such approach the number density of the QSOs 
at redshift z is the number density of corresponding 
peaks of density fluctuations (precursors of massive 
galaxies) which had collapsed between tc — tqso smd 
tc- These peaks we can see as QSOs at redshift z 
corresponding to time t(z) = tc{zc) + Ai: 



to 



tc{Zc) + At 



(4) 



The height of such peaks v is in the range of 
(vc, fc + Av), where 



Vc^l-69<yo'izc + l), 



Av^ 1.13(zc + l) 
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(5) 
(6) 
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as it follows from eq. (3). Then the number density 
of QSOs is 



nqsoiz) == a 



u^+Au 



p{v\vth,q) Npkiv) dv (7) 



where z, Vc and Az/ are defined by eq. 4, 5, 6 respec- 
tively. As we can see, the number density of QSOs 
iT-Qsoiz) depends upon the duration of the quasar 
stage tqso, the time delay At, parameters of the 
threshold function i^th and q, momenta of spectrum 
smoothed by galactic filter uo, cti, (T2 and fraction 
a of peaks-precursors passing quasar stage and visi- 
ble by Earth observer. If quasar stage is long-term 
{tqso ^ to) then the number density of QSOs at 
all z is equal to fraction of already collapsed peaks- 
precursors of massive galaxies visible from Earth at 
quasar stage of their evolution: ngsoiz) = ang{z). 
If quasar stage is short-term {tqso ^ tc), then it 
follows from (7) and (6) that 



nqsoiz) = l.lSa 



{Zc + 1) 
<^oto 



2.5 



Piv^^th^q) Npk{iy) TQSO- 

(8) 

It is valid for above mentioned estimation of tqso ~ 
10"^ - 10^ years. 

We suppose that the duration of the quasar stage 
TQSO J the time delay Ai and parameters of threshold 
function Uth and q do not depend on the z. Also we 
assume that the sample of bright quasars given by 
Schmidt et al. m and Boyle et al. p7| is complete for 
all z in the range 0.7-4.7. Here we do not discuss the 
problem concerned with the observational data but 
pay particular attention to the theoretical approach 
of their using for cosmological problems. 

Now we can calculate the number density of QSOs 
for set of redshifts and given values of the duration of 
QSO stage tqso, the time delay At, the fraction a 
of massive galaxies passing quasar stage, parameters 
Vth and q of threshold function and momenta of spec- 
trum ao, (7i, (72 • For present spectrum normalised in 
a certain way the momenta (Tq, ctij cr2 are calculated 
independently. The time delay At for given spectrum 
we find in following manner. Using eq.(8) we cal- 
culate the redshift of maximum Zmax of dependence 
iT-Qsoiz)- If Zmax is higher than one of observable 
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2.2 then we find the time delay 



Ai from equation 



driQsoiz) 



dz 



In those models where 



^obs 



2.2 we assume Ai 



= 0. (9) 

=2.2 

is less than or equal 
: 0. Rest four values 



cty Tqso, ^th and q are unknown. Finding of them 
by method of minimization of divergence vector is a 
matter of testing of known spectra. 



We try also to reconstruct the initial power spec- 
trum from observable data on number density of 
QSOs in wide range of redshifts and massive galaxies 
on z w 0. For this we consider momenta of spectrum 
(Joj fi: o'2 as unknown values and find them in the 
same way. 

In all cases we write cq.(7) for 6 points of data by 
Schmidt et al. on the zi = 4.7, Z2 = 4.05, Z3 = 3.7, 
Z4 = 3, Z5 = 2.8, zg = 2.2, for 3 points on curve by 
Boyle et al. Q on z-j — 2.0, zg — 1.5, zg = 0.7 and 
the eq. (2) for number density of massive galaxies on 
zio = 0: 



p(y; Vth, q) Npk{y) dv = n°^^„(z,), {i = 1, 2.. .9) 



p{v; Vth, q) Npk{v) dv = n°' 



bs 



(10) 

(11) 



where Vi, IS.Vi are calculated for Zi accordingly to 
(5) and (6) respectively. So, we solve the overdeter- 
mined system of equations (10), (11) (or (10) only) 
by method of minimization of the divergence vector 
using MathCad |22| software. 



3. Results 

3.1. Testing of spectra 

We shall test the preferable initial power spectra 
of density fluctuations on the flat Friedmanian back- 
ground expected in standard CDM model (^cdm = 
0.9, VLb = 0.1), tilted CDM ones and hybrid H+CDM 
one iytcDM ~ 0.6, flnDM = 0.3, fib = 0.1), which 
marginally match other observational data on large 
scale structure of the Universe [0-0. The trans- 
fer function T{k) in the spectra P{k) = Ak"T^{k), 
where A is constant of normalisation [J, is taken from 
the paper by Holtzman |2^. All spectra we nor- 
malised to produce the COBE data on CMB tem- 
perature anisotropy at angular scales 10° [gj, pSJ : 

< m' ' 

mula by Wilson and Silk |2^ for correlation func- 
tion of CMB temperature anisotropy measured by 
receiver with Gaussian response function and exact 
formula 1271, |2q] for correlation function C{a) =< 



>2= (1.1 ± 0.1) X 10"^ We calculate 
>= C(0; 10°) using the approximation for- 



^(0) ^(a) >. Both Sachs- Wolfe and Doppler ef- 
fects are taken into account (see also pa]). The con- 



^The parameter ti of transfer function fitting formula by Holtz- 
man psl are accepted equal to 1. 



tribution of tensor mode to AT/T is different in differ- 
ent models of inflation. In the most of them it is not 
dominating at the COBE angular scale (see review 
by |p9[ and references cited therein), in the some ones 
it is negligible even if n < 1 (for example "natural" 
inflation, [pO|). In this work, like ones cited here, we 
consider only models of inflation which give no sig- 
nificant tensor mode. But one can easily cross over 
to models of inflation with significant tensor mode 
when ratio ar = {AT/T)r^ / (AT/T)^ is known (here 
{AT/T)j, and {AT/T)g is contribution of tensor and 
scalar modes correspondingly). The constant of nor- 
malization then will be A' = A/{l+a'^), the momenta 
of spectrum a'^ = aj/ y'l + a^, peaks-precursors of 
QSOs will be collapsed at z^ = (zc + l)/-\/l + a"^ — 1, 
where values without (') are ones calculated with- 
out tensor mode. Parameters of cosmological models, 
constants of normalisation A of power spectra, their 
momenta ao, cti, 1T2 and biasing parameters b = l/ug, 
where as is the r.m.s. mass fluctuations in the top-hat 
sphere with radius 8h~^Mpc, are presented in Tabl.l 
. (For comparison, our normalization constant A for 
CDM model agrees with analogous constant from pa- 
per ll) reduced to h=0.5 and cr(10°) = 30.5/iX with 
accuracy better then 0.1%). 

For beginning we analyse the dependence of the 
ngsoiz) on tqso, At, vth and q for different spec- 
tra accepting a equal to 1. Results are presented in 
Fig.1-3.. 

As we can see, the number density of QSOs is pro- 
portional to TQSO, meanwhile the time delay displaces 
the maximum of uqso to lower z (Fig.l). Reducing 
of power spectrum at small scale results in steeper 
declination of nqsoiz) at high z and displacing of its 
maximum to lower z (Fig. 1-2). Increasing of the value 
of the threshold parameter q (tqso , i^th are constants) 
reduces the number density of QSOs at low redshift 
and displaces its maximum to high z (Fig. 3a). De- 
pendence nQso{z) on Vth {tqso and q are constants) 
is similar (Fig. 3b). Increasing of q with finding of 
corresponding i^th, which ensures Ug = n°'"*, results 
in steeper slope of ngsoiz) at low redshifts z < 2 
(Fig. 3c). The curve labeled by 5 in Fig. 3 is the 
number density of QSOs calculated without thresh- 
old function (vth = 0). For comparison we have cal- 
culated also the ngsoiz) without threshold function 
for another minimal linear overdensity corresponding 
to a bound object at z = 6^ = 1.33 adopted by 
Nusser and Silk ||lj| (fine labeled by 6 in Fig. 3a). 



Now let us find for each model such values of 
Oj Tqso, ^th and q which give minimal divergence 
iT-QSoiz) from observable one and ensure the num- 
ber density of massive galaxies at z « 0. In all cases 
analyzed here the amplitude of theoretical ngsoiz) 
is roughly equal to observable rig^Q (2.2) ~ 3 x 
10"^ h^Mpc~^ when tqso ~ 10'' years that is < to- It 
means that (8) is good approximation of (10) in the 
range 10"'^ < a < 1 and that ngso oc arqso- As 
the result, we can not find both values a and tqso 
simultaneously but their product arqso only. Ap- 
proximate solutions of system of equations (lO)-(ll) 
are obtained by using of MathCad software package 
and presented in the Table 2. The quasar abundances 
nqsoiz) for them are shown in the Fig. 4. 

As we can see any model does not explain the ob- 
servable nQgQiz) at high z. Therefore we have re- 
peated the same procedure without the equation for 
concentration of bright galaxies (11). Results are pre- 
sented in the Table 3 and Fig. 5. 

In this case the parameter a has other interpreta- 
tion: it is the fraction of peaks which are selected by 
threshold function piv^Vth t'P^^) and which pass 
through quasar stage and are visible by Earth ob- 
server. Let a' denote it. The ratio of the number 
density of all peaks selected in a such way to number 
density of bright galaxies 



K 






is calculated for each model and presented in Ta- 
ble 3. In all models, with the exception of standard 
CDM, K <^ 1 and can be interpreted as ai above 
mentioned, that is fraction of the peaks-precursors of 
bright galaxies passed through QSO stage. Thus, in 
these models the peaks resulting in bright galaxies 
through QSO stage are higher at average than main 
part of all peaks. The picture is inverse in standard 
CDM, where K ^ 7, and may mean that peaks of 
galactic mass AI > 2 x 10^^ h~^MQ survive after 
QSO stage and become bright galaxies if their gravi- 
tational potential holes are sufficiently deep. Whether 
it is valid or not can be proved only by means of elab- 
orate numerical simulations of the evolution of such 
peaks in different cosmological models. As it follows 
from Fig. 5, the theoretical number density of QSOs 
at low z is satisfactory in all models. At high z only 
CDM model with n = 0.7 gives the number density of 
QSOs sufficiently close to observable one. But even 
this model does not explain the observable number 



Table 1: The constant of normalisation of the spectra, theirs momenta cto, cti, 02 on galactic scale (i?/ = 
0.35/i^^il/pc) and biasing factor 6, defined by 6 = a^ , where erg is r.m.s. mass fluctuation on a scale of 8h^^ Mpc. 





CDM 


H + CDM 


n 


1 


0.8 


0.7 


1 


A 


7.33 X 10'^ 


1.62 X 10^ 


7.57 X 10^ 


7.51 X 10*' 


CTO 


4.78 


2.37 


1.67 


1.39 


0-1 


4.77 


2.18 


1.47 


0.91 


0-2 


7.93 


3.50 


2.31 


1.32 


h 


1.05 


1.77 


2.30 


1.66 
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Fig. 1. — Quasar abundance against redshift for different values of the tqso — 10^, 10^, 10^, 10®, 10^0 years (dashed 
lines 1, 2, 3, 4, 5 respectively in each panel) and time delay Ai = (a), Ai — 10^ years (b). At — 1.61 x 10^ years 
(c) for CDM spectrum with n = 1, ^cdm = 0.9 and Q,b = 0.1. The solid line with error bars is the comoving 
space density of quasars from Schmidt et al. (1991). The upper dashed line (5) is the number density of already 
collapsed peaks- precursors of bright galaxies. In all cases a = 1. 



density of QSOs at z > 4. Is it possible to explain 
the whole observable redshift distribution of quasar 
abundance given by Schmidt et al. M in such ap- 
proach, in principle? 

3.2. Finding of phenomenological spectrum 

For answer to this question we attempted to solve 
the system of equations (lO)-(ll) with respect to un- 
known values ctq, cti, CT2, cutqso, vtu and q. The time 
delay Ai is assumed to be equal to zero. Solutions 
are obtained in the same way and are presented in 
the Table 4. Quasar abundances ngsoiz) for them 
are shown in the Fig. 6. and match the observable 
data quite well. 



As we can see, the main feature of all solutions 
is the relation between momenta: cti < CT2 < fo- 
What kind of spectrum have such relations between 
its momenta on the galaxy scale ? We attempted to 
modify the standard CDM spectrum on small scale by 
force of reducing and enhancing of its power in order 
to obtain momenta of spectra from Table 4. Locations 
and amplitudes of such modifications were unknown 
parameters and were found from 3 equations for given 
momenta ci, 172, fo using MathCad software package. 
It was successful for momenta from the column 1 in 
Table 4 only, for which parameter 7 = (T^/crocr2 = 
0.3. The phenomenological power spectrum of density 
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Fig. 2. — Quasar abundance against redshift for CDM spectrum with flcDM = 0.9, ilf, = 0.1, n ~ 0.8 (a), n — 0.7 
(b) and hybrid H+CDM model {nnDM = 0.3, ncoM = 0.6, fit = 0.1) for different values of the tqso = 10^, 10^ 
10^, 10*, 10^° years (dashed lines 1, 2, 3, 4, 5 respectively in each panel). Time delay At is 5.2 x 10* years (a), 
(b) and (c). The rest is the same as in Fig. 1. 

Table 2: Solutions argso, t^th and q of system of equations (lO)-(ll) for different spectra. 





CDM 


H + CDM 


n 


1 


0.8 


0.7 


1 


aTQsoiyears) 
At (years) 


3.54 X 10^ 

1.55 X 10^ 
5.22 
1.95 


4.56 X lO'^ 
5.19 X 10* 

2.73 

6.60 


2.39 X 10^ 
0.00 
2.52 
12.4 


3.82 X lO'' 
0.00 
2.37 
4.30 



fluctuations obtained in such way is 



Pph(k) — AkTfjjjj^j 



1, 

gl-Cfc/fco)"^ 



k < kg, 

ko < k < ki, 
fc > fci, 



where ko = 0A8hMpc-\ ki = 0.82hMpc-^, p = 2.3, 
A is the same as in standard CDM. The spectrum 
has shelf-like reducing of power at /c ^ 0.7hMpc~^ 
in comparison with standard CDM normalised to the 
COBE 10° angular scale temperature anisotropy of 
CMB {Piph{k) spectrum in Fig. 7). 

The correlation function of galaxies for it does not 
contradict observable data. Other cosmological con- 
sequences of such spectrum are being investigated yet. 

The parameter 7 for momenta from columns 2-5 in 
Table 4 is extremely small (0.04-0.0004), so that phe- 
nomenological power spectrum for them, obviously. 



must be extraordinary too. Indeed, spectrum with 
momenta from column 2 (P2 ph{k) spectrum in Fig. 7) 
has two narrow bumps at fc « 0.02 and « lOh Mpc^^ 
and shelf-like reducing of power at fc ~ 0.03hMpc~^ 
in comparison with standard CDM. Since its corre- 
lation function of bright galaxies strongly contradicts 
observable one, the solutions from columns 2-5 in Ta- 
ble 4 can be interpreted as meaningless from physical 
point of view. So, only spectra with momenta from 
column 1 (7 « 0.3), or close to them, are cosmologi- 
cally significant. 

3.3. Estimation of fraction of host galaxies 

Masses of galaxies which pass through quasar stage 
(host galaxies) are Mg > 2 x IO^^H'^Mq, so that it 
is coordinated with the black hole accretion models 
of quasar phenomenon flS, Efl, R^, R3| . An optical lu- 
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Fig. 3. — Quasar abundance against redshift in tilted CDM model {n — 0.7) for different parameters of threshold 
function: (a) - g = 16 (line 1), q = 12 (2), q = 8 (3), q — 4 (4), without threshold function (5), tqso — 6.9 x lO"^ 
years, i^th — 5.89; (b) - vth — 10 (line 1), i^th = 8 (2), i^th = 5.89 (3), i^th = 3 (4), without threshold function (5), 
TQSO = 6.9 X 10^ years, q = 8.08; (c) - g = 16, m. = 2.5 (line 1), q = 12.4, i^th = 2.52 (2), 9 = 8, vth = 2.6 (3), 



g = 4, vth = 2.95 (4), without threshold function (5), argso 
eq. 11). The rest is the same as in Fig. 1. 



2.4 X 10 years, {lyth for given q was obtained from 



minosity of the bright quasars is L > 10''^ h^^ erg/ s. 
Their lifetime multiplied by a in all cases is a tqso ^ 
10^ years, because we denote tqso = 10^ T^/a, where 
r5 is a TQSO in the units of 10^. If the fraction of 
the black hole mass converted into optical radiation 
is e, the ratio of the quasar's mass to that of the host 
galaxy is F, then 



Afg w 2 X 10^ L4 



TrJia eF) h- 



that gives 



a « 10 ■^ L47 T-5/(eo.i Fo.oi) h '' 

where eo.i = e/0.1, Fq.oi = F/0.01, L47 = L/10''^. 

Thus, in the models with phenomenological power 
spectrum and spectra tested here only less than 1% 
of massive galaxies pass through quasar stage and are 
visible by Earth observer. 

4. Discussion 

Recently we |3J, Q carried out the similar calcu- 
lations of evolution of quasar number density for the 
same spectra (CDM with n=l, 0.8, 0.7 and H+CDM) 
normalised to ag = b~^, where galactic biasing pa- 
rameter bg was calculated in the framework of Gaus- 



sian statistics of density peaks EI 



and was equal 



to 1.40, 1.47, 1.53 and 1.76 respectively. From com- 
parison of results of this work with previous ones it 
follows that the dependence of quasar number den- 
sity on z is very sensitive to the normalisation of 
spectrum in all models. In our previous papers the 
modified on small scale standard CDM spectrum for 
reproducing of the observable data on number den- 
sity of QSOs at different redshifts has reducing of 
power on scale 1 < A: < 10hMpc~^ and bump on scale 
k ^ 10hMpc~^ (dotted line in Fig. 7). Its difference 
from spectrum shown in Fig. 7 is caused by different 
power of initial CDM spectrum on k < lhMpc~^ as 
a result of different normalisation. 

It should be remarked here, though formally this 
method is sensitive to spectrum on small scale up to 
k ~ lOh Mpc^^ when it is smoothed by Gaussian 
filter function with Rf ~ 0.35/i^^A/pc (or top-hat 
Rth — 0-6h~^Mpc) nevertheless we are not sure that 
calculation of peak number density of galaxy scale is 
exact for cases of spectrum like standard CDM and, 
especially, CDM + short wave bump models. This 
is concerned mainly with the well known problem 
'cloud-in-cloud' unresolved by now in the frame of 
peak formalizm used here. 

Now let us compare our results with those by 
Nusser and Silk @. The curve nQso{z) for 6° = 1.33 



Table 3: Solutions a tqso^ ^th, Q of system of equations (10) and ratio K of number density of peaks which has 
passed through quasar stage to one of bright galaxies. 





CDM 


H + CDM 


n 


1 


0.8 


0.7 


1 


a' TQSO (years) 


5.73 X 10^ 


1.83 X 10^ 


6.91 X 10^ 


1.64 X 10^ 


At (years) 


1.3 X 10^ 


5.19 X 10« 


0.00 


0.00 


t'i/i 


1.09 


7.26 


5.89 


4.95 


q 


4.96 


4.85 


8.08 


12.32 


K 


7.2 


0.02 


0.006 


0.003 



Table 4: Partial solutions aTqso, i^th, q, cq, cti and CT2 of system of equations (lO)-(ll) 



aTQSO (years) 


3.20 X 10^ 


3.51 X 10^ 


3.53 X 10^ 


4.41 X 10^ 


4.0 X 10* 


At (years) 


0.00 


0.00 


0.00 


0.00 


0.00 


I'th 


2.06 


1.69 


1.69 


1.67 


1.67 


q 


8.71 


9.28 


9.28 


9.31 


9.31 


o-Q 


2.38 


2.77 


2.78 


2.81 


2.82 


0-1 


1.48 


0.22 


0.17 


0.02 


0.002 


0-2 


2.37 


0.43 


0.34 


0.03 


0.004 



(line 6 in Fig. 3) without threshold function, recalcu- 
lated to aTQso — 10* years, coordinates with the cor- 
responding curve in Fig. 2 of work by Nusser & Silk 
flj] quite satisfactory at high z. The some divergence 
of these curves at low z are caused mainly by diferent 
ways of definition the mass of peaks. 

In our work the value d° — 1.68 corresponding to 
collapse of spherical-symmetrical peaks, is adopted. 
Therefore, in our approach the QSOs are formed in 
spheroidal peaks only. The fraction of spheroidal 
peaks is < 4% [Q, the fraction of massive galax- 
ies with M > 2 X 10^^ h~^MQ passed through quasar 
stage in all cases does not exceed 1%. So, the as- 
sumption that QSOs appear only in spherical peaks 
of density fluctuations is not contradictory in our ap- 
proach. But its conformity to reality should be proved 
or objected by elaborated numerical simulations of 
the evolution of QSOs which start from initial condi- 
tions determined at the linear stage of the evolution 
of peaks. Substantiations of the value S° as well as 
dependence or independence of both quasar lifetime 
TQso and time delay At on redshift z will be the sub- 
ject of the future works, c We put constraints on 



the cosmological models based on their ability to re- 
produce the rise of quasar abundance up to redshifts 
w 2.2 and following fall up to « 4.7 assuming that 
quasar nature is the same at different redshifts and 
TQSO J At and a do not depend on z. Moreover, if 
we have only 3 unknown values aTqso, i^th and q in 
equation system (10-11) then we can not reproduce 
the nQ%Q{z) given by Boyle et al. [1^] and Schmidt 
et al. iy with sufficient accuracy in the models tested 
here. But when cro, fi and (T3 are supposed to be un- 
known values too, such reproducing is possible. Ap- 
parently, for other models of quasars, which explain 
the rise and fall of the nqgQ^z) by astrophysics of 
quasars (e.g. [0), the constraints on the cosmolog- 
ical models may be other. But it is the subject of 
separate work. Here we only confirm that such de- 
pendence of quasar abundance on z can be explained 
in the framework of most simple astrophysical models 
of quasars, which look rather plausible and motivated 
in order to be put in the base of similar research. 

The phenomenological power spectrum obtained 
here (§3.2) has shelf-like reducing of power at fc ~ 
0.7h Mpc~^. Such spectra with broken scale invari- 
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Fig. 4. — Quasar abundance against redshift ob- 
tained by minimization of divergence vector including 
also number density of bright galaxies at z ~ 0: 1 - 
CDM n = 1, 2 - CDM n = 0.8, 3 - CDM n = 0.7, 4 - 
H+CDM n = 1. 



ance is generated in some inflationary scenarios (e.g. 
m, H, 0, |ll). The typical scale and the height of 
the shelf-like reducing determined here can be con- 
nected with some constants, which characterise the 
underlying inflationary models. 

If we treat the used observational data more crit- 
ically and suppose that nQsoiz) given by Schmidt 
et al. M and Boyle et al. jl^] is only lower limit 
for number density of QSOs at all z then the next 
constraints follow from the presented results: on the 
small scale the initial power spectrum of density fluc- 
tuations is nearly CDM with n > 0.7, bg < 2.3 and 
aTQso > 5 X 10^ years, Ai < 5 x 10* years. 

5. Conclusions 

Thus, the main assumption that quasar phenomenon 
is active short term stage of evolution of some small 
fraction of massive galaxies with M > 2x10^^ H^^Mq 
and appears in peaks of random Gaussian density 
fluctuation field allows to explain the general feature 
of riQ^gQiz) given by Boyle et al. [0 and Schmidt et 
al. W: fast increasing of number density from z ss 
to z ~ 2.2 up to value ~ 3 x 10""^ h^Mpc'^ and 
following monotonous decreasing on z > 2.5. Appar- 
ently, a constant quasar abundance from z « 2 to 4 



Fig. 5. — Quasar abundance against redshift ob- 
tained by minimization of divergence vector without 
abundance of bright galaxies: 1 - CDM n = 1, 2 - 
CDM n = 0.8, 3 - CDM n = 0.7. 4 - H+CDM n = 1. 



|(42| can not be explained under these assumptions. 

If the nature of bright QSOs is the same at dif- 
ferent redshifts {tqsoi ^^i oi are constants) then any 
model spectrum tested here such as standard CDM 
spectrum [ilcDM = 0.9, Q.b — 0.1), tilted one with 
n = 0.7,0.8 and hybrid H+CDM {VLcdm = 0.6, 
^HDM = 0.3, Qb — 0.1, n — I) explains the number 
density of such QSOs in redshift range z « 0.5 — 2.2, 
but in range z « 0.5 — 3.5 only tilted CDM with 
n = 0.7 {tqso « 7 X 10^ years, ai « 0.01, At = 0) 
does explain, and in range z « 0.5 — 4.7 - none of 
them. Consequently more crucial for testing of spec- 
tra is number density of QSOs on high z. 

Exact reproducing of number density of QSOs 
a,t z — 0.5 — 4.7 under above mentioned assump- 
tions can be done in the flat cosmological model 
with the density fluctuation power spectrum which 
gives (To ~ 2.4 — 2.8 and 7 < 0.3 at galaxy scale 
{Rf = 0.35 /i-i Mpc). It predicts the QSO's lifetime 
'Tqso ~ (3-5 — 4.5) x 10^/a years, where a is the frac- 
tion of massive galaxies passing through quasar stage. 
If physically motivated quasar lifetime is ~ 10^ — 10® 
years, then a w 0.01 — 0.001 that matches the esti- 
mations of a done above. 

Thus, redshift distribution of quasar abundance is 
very sensitive to the amplitude and slope of the initial 
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Fig. 6. — Quasar abundance against redshift for 
found momenta of spectra, parameters of threshold 
function and argso- 

power spectrum of density perturbations on galaxy 
scale. But it can be effective test when complete num- 
ber density of QSOs with mass larger than fixed value 
will be known with confidence for all redshifts and 
physical parameters connected with nature of quasar 
phenomenon (the quasar lifetime tqso, time delay 
At, fraction of galaxies which pass through quasar 
stage a, linear amplitude of peaks 6° collapsing just 
now and threshold function) will be substantiated. 
The last needs the detailed investigation of connec- 
tion of quasar phenomenon with initial configuration 
of density fluctuation by the help of numerical simu- 
lations of its evolution. 

We acknowledge useful discussions with B. Hnatyk, 
V.N. Lukash and B.V. Romberg. 
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